Three dimensional complex plasma structures in a combined radio frequency and 

direct current discharge 
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We report on the first detailed analysis of large three dimensional (3D) complex plasma structures 
in experiments performed in pure rf and combined rf+dc discharge modes. Inductively coupled 
plasma (ICP) is generated by an rf coil wrapped around the vertically positioned cylindrical glass 
tube at a pressure of 0.3 mbar. In addition, dc plasma can be generated by applying voltage to 
the electrodes at the ends of the tube far from the rf coil. The injected monodisperse particles 
are levitated in the plasma below the coil. A scanning laser sheet and a high resolution camera 
are used to determine the 3D positions of about 10 5 particles. The observed bowl-shaped particle 
clouds reveal coexistence of various structures, including well-distinguished solid-like, less ordered 
liquid-like, and pronounced string-like phases. New criteria to identify string-like structures are 
proposed. 
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I. INTRODUCTION 

Complex (dusty) plasmas - systems consisting of 
highly charged micron-size particles in a neutralizing 
plasma background - exhibit an extremely rich variety 
of interesting phenomena |l|, [2j, including phase tran- 
sitions 0-0] ' solid- liquid interfaces and crystallization 
front propa gati ons |j|, fluid-like instabilities Q, noctilu- 
cent clouds [lfl], etc. Amongst these formation and evo- 
lution of three-dimensional well-ordered structures are 
of particular interest p], llTH18| . This is mainly because 
high temporal and spatial resolution allows us to inves- 
tigate these structures along with various related phe- 
nomena at the individual particle level. Fully resolved 
quasi-atomistic (undamped) particle dynamics provides 
new insight into natural atomic and molecular systems, 
whose dynamics cannot be resolved in such detail. 

Since the discovery of plasma crystals [3j-|6| many ex- 
periments have been performed to quantify the differ- 
ent states of complex plasmas by analyzing the scattered 
laser-light of the microparticles that were introduced in 
the plasma. Most of these experiments were related to 
two dimensional (2D) systems (e.g. pUa 181 fl9|). Only 
a few experiments were focused on the analysis of large 
(~ 10 5 particles) 3D dusty structures @, [TlTlrl | mostly 
performed in parallel plate rf discharges. 

Here we report on the first fully 3D reconstruction 
of a large microparticle cloud in inductively coupled 
electrodeless rf (or rf+dc) plasma in cylindrical geometry 
and detailed analysis of 3D structural data. In contrast 
to other typical complex plasma experiments in dc 
discharge tubes (see e.g. [20]) the particles in the present 
experiment are not levitated in standing striations, 
where strong variations in the electric field (and other 
related plasma parameters) lead to inhomogeneities 
in particle clouds even on small scales. Instead, the 
particles are confined in a vertically mounted discharge 
tube, where levitation becomes possible due to an 
rf (or rf+dc) electric field and an upward neutral gas 
flow balancing the particle gravity and the ion drag force. 
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FIG. 1: (Color online) Sketch of the setup and of the force 
balance on the particles levitating in the combined rf and dc 
plasmas where Fi is ion drag force, F g gravity, Fdc electric 
force by dc plasma and Fqf is neutral drag force (due to the 
gas flow). 



II. EXPERIMENT 

The experiment is performed in the "PK-4" setup de- 
signed as a laboratory prototype of the next generation 
complex plasma experimental facility onboard the ISS. 
The heart of this setup is a vertically oriented "U"- 
shaped glass tube of 30 mm in diameter and 300 mm 
working length. The setup is equipped with all stan- 
dard systems for operating the plasma, vacuum pumps, 
flow controllers, high voltage supplies for dc and rf driven 
plasmas, and a system for particle observation: illumina- 
tion laser spread into a thin sheath and a fast camera 
perpendicular to it [21] . Specific for this experiment is 
the system for particle scanning, where illumination laser 




FIG. 2: (Color online) Dust structures observed in the pure rf 
(left panel) and combined rf+lmA dc (right panel) discharge 
with general (a) and side (b) views. Particles are color-coded 
by their radial distance from the tube axis (the distance is in 
mm units). The white squares show relatively homogenous 
regions used for the structural analysis. The lines indicate 
approximately the local orientations of the (string- like) struc- 
tures. Considered structures consist of about 10 5 particles, 
which is comparable to the biggest complex plasma systems 
analyzed so far. 



and camera are moved together by a motor for scanning 
through the depth of the particle cloud, described previ- 
ously in [lj] . Figure 1 shows the sketch of the experiment 
indicating the position of the particle cloud and the main 
forces acting on the particles. 

The inductively coupled plasma is sustained by an rf 
current at frequency of 81.36 MHz driven through the 
wire looped around the middle of the tube. The in- 
put power of the rf coil is 1.6 W, creating an elliptically 
shaped plasma. In addition, plasma can be generated 
through the whole tube by applying a dc voltage to the 
electrodes at its ends. Experiments are performed in a 
mixture of argon (90%) and oxygen (10%) at a pressure 
of 0.3 mbar and an upwards gas flow of 0.25 seem. Al- 
though similar global behavior is observed in pure argon 
and neon plasmas full structural analysis was done for 
this Ar:0 gas mixture. Gas flow is used in order to avoid 
the plasma degradation and to provide additional com- 
pensation for the particle gravity. This additional com- 
pensation is important to position the particle cloud in 
the more homogeneous plasma, where the particle system 
is less disturbed by plasma-related instabilities. 

Melamine formaldehyde microparticles of 2.55 /im in 
diameter are injected from the top of the tube. Due to 
the balance of the upward electrical and neutral drag 
(associated with the gas flow) forces with the downward 



gravity and ion drag forces, particles can levitate in a 
certain region below the rf coil (see Figs. Q] and H). 

In the basic experiment particles are levitated in pure 
rf plasma with the mentioned constant gas flow. Ap- 
proximately a minute after the particle injection a scan 
through the particle cloud is performed with 26 frames 
per mm (and 33 frames per second). After the scan a 
dc voltage is applied to the electrodes, creating a cur- 
rent limited plasma. The dc electric field is oriented to 
additionally compensate the particle weight. The scan 
through the particle cloud is repeated for the new plasma 
conditions. Current is changed to 0.2, 0.3, 0.4, 0.6, 1 
mA and particles positions are determined for each cur- 
rent. The third dimension (perpendicular to the field of 
view) of particles positions are reconstructed based on 
the intensity profile of the illumination laser as described 
in [14|. During scans every particle position is recon- 
structed in 7 to 10 frames from which we determined the 
average particle position in all three dimensions. The 
average vibrational motion of particles during its detec- 
tion is in most cases limited to one or two pixels, cor- 
responding to 0.02 to 0.04 mm which is much smaller 
compared to the average interpaticle distance of around 
0.2 mm. These numbers for particle displacements are in 
good agreement with our estimates of particle diffusion 
for parameters relevant to the present experiment. The 
Epstein damping rate 22j, characterizing frictional dis- 
sipation in grain-neutral collisions is i/f r ~ 130 s _1 . The 
dimensionless parameter, which can appropriately char- 
acterize the damping strength is the so called "damping 
index" [23j . £ = U{ T Ay^M/T, where A is the interparti- 
cle distance, M is the particle mass, and T is the kinetic 
temperature of the particle system [in the following we 
assume that T is equal to the neutral gas (room) tem- 
perature]. Physically, £ is the ratio between the charac- 
teristic interparticle spacing and mean ballistic free path 
of the particle. The particle motion is merely ballistic 
(Newtonian) for £ <§C 1, while for £ ^ 1 it is overdamped 
(Brownian). For our parameters we get £ ~ 40, which in- 
dicates that the particle motion approaches overdamped 
regime. The characteristic value of the diffusion coef- 
ficient for the system of particles containing solid and 
fluid phases can be roughly taken as the diffusion co- 
efficient at the fluid-solid phase transition, D ~ O.lDo, 
where Do = T/(Mvf T ) is the bare Brownian diffusion co- 
efficient (for non- interacting particles). This relation is 
known as the dynamical criterion for freezing |24| . In our 
conditions Do~3x 10~ 5 cm 2 /s and hence D~3x 10~ 6 
cm 2 /s. The characteristic displacement during r ~ 0.3 s 
(10 frames) is £ ~ \JQDt ~ 0.02 mm, in good agreement 
with the observations. The characteristic time required 
to travel one interparticle distance, ta ~ A 2 /6D ~ 25 s, 
is longer than the scanning time. 

Dust structures observed in pure rf plasma have a 
bowl-like shape. With the applied dc voltage the par- 
ticle density near the tube center decreases, and a void 
is gradually formed when the dc current increases (see 
Figure [5]and|3]). We propose the following explanation 





FIG. 3: (Color online) Axial crossections of the particle cloud 
levitated in combined rf and dc discharge with different dc 
currents (a) 0.1 mA, (b) 0.2 mA, (c) 0.4 mA and (d) 0.6 mA. 
It is clearly visible that the global distribution of the particles 
is systematically effected by the dc current increase, slowly de- 
creasing the number of particles levitating along the center of 
the tube and gradually forming the void around them. Par- 
ticles are color coded by the radial position indicated by the 
color (distances are in mm units). 



for this observation. The application of dc voltage shifts 
the maximum of the rf plasma production down with re- 
spect to the position of the rf coil (i.e. in the direction 
of the dc electric field) . This is evidenced by the analysis 
of the distribution of the integral plasma glow along the 
tube axis, obtained using a video imaging. This implies 
that the plasma density increases below the coil (and de- 
creases above it) in comparison to the pure rf generation 
case. The result is that the relative importance of the 
ion drag force in this region increases. Increasing the ab- 
solute magnitude of the ratio of the ion drag force to the 
electric force (which should not necessarily exceed unity) 
has a consequence that the particle levitation becomes 
impossible in some vicinity around the tube axis. Very 
close to the tube axis some particles can still levitate due 
to the upward neutral gas flow, which has a parabolic 
profile and produces maximal effect on the axis. Levita- 
tion remains also possible sufficiently far from the axis, 
where the plasma density decreases and the total ion flow 
velocity increases due to the additional radial component 
of the discharge electric field (both effects lower the rela- 



tive importance of the ion drag force 25]). For interme- 
diate radial position no levitation is possible and a void 
is formed. This explanation is in qualitative agreement 
with observations, cf. Fig. [5]and|31 



FIG. 4: (Color online) The typical particle distribution on the 
plane of bond order parameters q^ and qe, (for pure rf plasma). 
Particles are color-coded by the wa value. Rotational invari- 
ants q±, q§, u>4 were calculated by using 12 nearest neighbors; 
those for perfect fcc/hcp/ico lattice types are also indicated. 
The distribution reveals presence of solid-like (hep-like and 
fee- like), weakly disordered (glassy), liquid- like and string- like 
phases in the cloud. Inset shows the radial distribution func- 
tion (RDF) for rf (solid red line) and rf+lmA dc (solid green 
line) plasma. The cumulative distributions N(< r) (mean 
number of particles located inside sphere of radius r in mm 
units) of the RDFs are shown by dashed lines. Both the cumu- 
lative distribution N(< r) and splitting of the first maximum 
of the RDF indicate the presence of string-like structures. 



III. STRUCTURAL PROPERTIES 

Irrespective of the applied dc voltage, the observed 
clouds of particles are not very homogeneous. Even in 
a pure rf regime, typical interparticle separations in the 
central part of the cloud can exceed those in the periph- 
eral regions close to the cloud boundaries by a factor of 
about two. For this reason, a relatively small peripheral 
part of the cloud sketched in Fig. Hfb) has been chosen 
for the detailed analysis of the particle systems structural 
properties. This part is sufficiently small (especially its 
radial extent) so that the system inside is reasonably ho- 
mogeneous. At the same time, it contains enough parti- 
cles (~ 10 4 ) to yield reasonable statistics. For the rest of 
the article we compare the structural arrangement of the 
levitated particle cloud under two extreme conditions, 
pure rf plasma and combined rf and dc plasma with max- 
imal dc current of 1mA. 

To determine the local structural properties of the 
three-dimensional particle system we use the bond or- 



der parameter method (261 127] . which has been widely 



used to characterize order in simple fluids, solids and 



glasses, hard-sphere systems, colloidal suspensions, 3D 
complex plasmas (see e.g. Refs. [121 EM [28| and refer- 
ences therein). In this method, the rotational invariants 
of rank I of both second qi{i) and third wi(i) order are 
calculated for each particle i in the system from the vec- 
tors (bonds) connecting its center with the centers of the 



N nn (i) nearest neighboring particles: 
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where qi m (i) = N an (i) 1 J2f=i Y lm(rij), Y im arc the 
spherical harmonics and r^ = r, — Yj are vectors connect- 
ing centers of particles i and j. In Eq.© 

,— ' mi rri2 m3 

denote the Wigner 3j-symbols, and the summation in 
the latter expression is performed over all the indexes 
mi = —I, ..., I satisfying the condition mi +m2 + m^ = 0. 
The calculated rotational invariants q,. wj. are then com- 
pared with those for ideal lattices [7], [2|| (27|. Here, we 
are specifically interested in identifying face-centered cu- 
bic (fee), hexagonal close-packed (hep), icosahedral (ico) 
and body-centered cubic (bee) lattice types and, there- 
fore, use the invariants q±, qe, IU4, wq calculated using the 
fixed numbers of N nn = 12 (fcc/hcp/ico) and N nn = 8 
(bec) nearest neighbors, respectively. A particle whose 
coordinates in the 4-dimensional space ((74, qe, W4, u>e) are 
sufficiently close to those of the ideal fee (hep, ico, bec) 
lattice is counted as fcc-like (hep, ico, bcc-like) particle. 

An example from this analysis is presented in Fig. 21 
which shows representative particle distributions on the 
plane ((74, q§ ) taken for the particle cloud observed in 
pure rf plasma. The distribution of particles versus rota- 
tional invariants 94, qe reveals the presence of solid- like 
(mostly hep-like), string-like (big values of q^ and 104), 
medium ordered glassy (0.4 < g 6 < 0.44) and weakly or- 
dered liquid structures (54 < 0.4 and qe < 0.4 ) in the 
analyzed part of the cloud. The corresponding domains 
are indicated in the plot. 

The inset in Fig. U shows the radial distribution func- 
tion (RDF) g(r) for both pure rf (solid red line) and com- 
bined rf+lmA dc (solid green line) plasma. Cumulative 
distributions iV(< r) (mean number of neighboring par- 
ticles in the sphere of radius r) for these RDFs are also 
plotted (dashed lines). Both g(r) and N(< r) indicate 
the presence of the string-like structures (SLSs) in the 
particle cloud: The splitting of the first RDF maximum 
(at N ~ 12) is one of the indicators of developed SLSs in 
the system (see |15l|). 

In view of the ubiquitous character of string-like struc- 
tures in various complex plasma experiments (e.g. jl|, 




FIG. 5: (Color online) Typical particle distribution on the 
plane of the order parameters 54 — qe (calculated by using 6 
nearest neighbors) . The domains occupied by both the quasi- 
isotropic and string-like structures (consisting of JVsls > 7 
microparticles) are indicated. The color coding from blue 
to red corresponds to the increasing probability of finding a 
particle in a given region of (54, qa) plane. 



|29|). including the present one, it would be of value to de- 
velop simple and convincing indicators that can be used 
to identify the transition from isotropic to anisotropic 
(SL) structures. Here we put forward and briefly discuss 
three related possibilities. 

First, we consider a new descriptor of the SL struc- 
tures, which is related to the distribution of particles in 
the plane (54, q§), where the rotational invariants (74, qe 
are calculated using a fixed number of 6 nearest neigh- 
bors (instead of 12 and 8 to identify isotropic crystalline 
structures considered above). This choice of the num- 
ber of nearest neighbors used for detection of string like 
structures is to some extend arbitrary. In the analyzed 
case, observation of 6 nearest particles yields good statis- 
tic and clear identification of "good" strings: In an ideal 
case when all 7 particles tend to lie on the same line, 
the value of 54 tends to unity. Using smaller number of 
nearest neighbors results in higher values of 54 and qg. 
Nevertheless, the pronounced peak in the 94 values, well 
separated from the values characterizing the rest of the 
system would give clear identification of the anisotropy. 
A typical example of such distribution (for the regime of 
pure rf plasma) is shown in Fig. [5) Clearly, the existence 
of SL structures (consisting of -/Vsls ^ 7 microparticles) 
can be easily identified. 

Another useful measure of the system anisotropy is 
the distribution function of the distance r m ; n between a 
test particle and its nearest neighbor. The corresponding 
probability distribution functions (PDFs) are shown in 




FIG. 6: (Color online) Observed string-like structures in pure 
rf (a) and rf+lmA dc (b) plasma. Particles are color-coded by 
the parameter r m i n (mm) which is the distance to the nearest 
neighbor. Inset shows corresponding probability distribution 
functions. The cumulative distribution (green line) are also 
plotted. The PDF reveal shoulder-like behavior at distances 
r far less than the averaged interparticle distance, indicating 
the presence of string-like structures in the system. 




FIG. 7: (Color online) Probability distribution function P(a) 
in an rf (solid red line) and rf+ rf+lmA dc (solid green line) 
plasmas. Similar PDFs are alsoplotted for isotropic liquid- 
like (the data set is taken from [J, solid blue line) and solid- 
like (the data set is taken from [la |. solid orange line) 3D 
complex plasmas. Corresponding cumulative distributions are 
also plotted by dashed lines of the same color. The inset shows 
the distribution of the string-like structures over the number 
of particles they contain for rf (solid red line) and rf+lmA dc 
(solid green line) complex plasmas. 



tions P{a) for the observed complex plasmas for rf (red 
line) and rf+lmA dc (green line) plasma. Additionally, 
we plotted P(a) for the solid-like and isotropic liquid- 
like complex plasma systems (calculated from the data 
of Refs. Q and [HI]) to emphasize structural anisotropy 
observed in the present experiments. The inset in Fig. [7J 
shows distributions of the observed SLSs over the num- 
ber of particles they contain. These distributions clearly 
reveal a power law-like behavior, indicating the absence 
of any spatial scales characterizing the anisotropy of the 
particle systems under consideration. 



Fig. [5] together with the overview of the particle clouds. 
The presence of the SL structures is evidenced by a non- 
gaussian behavior of these PDFs. Thus, the shape of 
the PDFs of r m j n turns out to be a useful tool to identify 
system anisotropy, in addition to the widely used scaling- 
index method [lji [3(| . 

An important characteristic of SL structures is the dis- 
tribution over the number of particles they contain. To 
find this distribution we define two nearest neighbors for 
each particle i, their positions are characterized by r^i, 
r^2 and r,, respectively. We then calculate the angle a 
between the vectors r^ — r,-,i and r,,2 — r^. If a is suffi- 
ciently small (a < a*), it is postulated that the particles 
belong to the same string (in practice a* sw 10° is used, 
the sensitivity to the exact value of a* is rather weak). 
Clustering algorithm (see Q) is used to define the to- 
tal number of particles in the string and to obtain the 
related distribution of SL structures over the number of 
particles. Figure [7J shows typical probability distribu- 



IV. CONCLUSION 



To conclude, we present first detailed analysis of large 
3D complex plasma structures in rf and combined rf+dc 
generated plasmas. Structural analysis reveals a rather 
complicated mixture of different phases in the observed 
complex plasmas: presence of solid-like, liquid-like and 
string-like structures in the system is well established. 
Variations of the dc current significantly change the 
global structure of the particle cloud, but do not affect 
the local structural properties significantly. Finally, we 
propose new indicators to identify and quantify string- 
like structures appearing in complex plasmas. One of 
the observations from the present experiment is that the 
identified SLSs exhibit scale-free power-law-like distribu- 
tions over the number of particles they contain. 
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